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1. INTRODUCTION 

The drug delivery process of general anaesthesia is most commonly described by the 
Pharmacokinetic/Pharmacodynamic (PK/PD) model. The PK model is a linear time invariant (LTI) model 
that relates the drug dosage to drug concentration in blood plasma, while PD model is a static nonlinear 
model that relates the drug concentration in the blood plasma to the resulting drug effect [1]. 

One important property of the PK model is that it is a positive system; the state variables (drug 
concentration in each compartment) can never be less than zero. Imposing positiveness in the design of 
closed-loop system allows one to simplify the stability analysis [2] and reduces the problem to a real bounded 
uncertainty gain problem [3]. Therefore, the control of general anaesthesia can be treated as a stabilization 
problem of a positive system. 

Closed-loop stabilization of positive linear system has been considered in several previous works. 
For example, Kaczorek [4] proposed a sufficient condition for state feedback controller based on Gersgorin's 
theorem and quadratic programming. Later, Leenheer and Aeyels [5] solve the stabilisation of a state 
feedback using algebraic approach. Boyd et al [6] and Gao et al [7] proposed a necessary and sufficient 
condition to stabilised the closed-loop positive linear system formulated as Linear Matrix Inequality (LMI). 
Rami et al [8-9] also proposed a feedback stabilisation based on linear programming (LP) method. 

The LP approach offers several advantages over the other methods: it provides both the necessary 
and sufficient condition to guarantee asymptotic stability of the system, it has a lower computational 
complexity compared to quadratic programming or LMI, and it can be extended easily to uncertain plants. 
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Moreover, the LP inequalities can be used to design a positive interval observer [10]. In anaesthesia 
control, the drug concentration in each compartment was not possible to be measured in real time. Hence, an 
observer is needed to estimate the states. However, the classical Luenberger observer does not guarantee the 
positiveness of the estimated states. This can be achieved by using the LP approach that guarantees the 
estimated states to remain in the nonnegative orthant. 

Conventionally, state feedback controllers have been designed using the LMI [11] or pole placement 
[12-13] methods. However, for problem with only partial information on the state is available, static output 
stabilization problem has not been formulated and solved exactly as an LMI problem or pole placement 
approach [8]. Furthermore, it can be difficult to use pole placement method for positive system since it is not 
even know which poles are desireable to ensure positivity of the closed-loop system. Hence, the LP approach 
has the advantage of easy handling of static output stabilization problem for positive systems. 

In this paper, uncertainty was assumed to be occurred in PK model. LP approach was used to design 
the state-feedback controller and the positive interval observer, by taking into account the possible variability 
of parameters in the PK model. 


2. RESEARCH METHOD 

Several important notations are stated here. R"; denotes the nonnegative orthant of the n- 
dimensional real space R®. M™ denotes the transpose of the matrix M. Each elements of matrix M is 
represented by mj, where i is the number of row and j is the number of column. Mj; > 0 means that all the 


elements are nonnegative. M and M represents the upper bound and lower bound of the uncertain matrix M. 


B= Ax+ Bu 
y=Cx (1) 


a. Plant model 


The PK model can be described as a four-compartment model and expressed in the following state 
space Equation: 


[ V. V. 
-k,-k,,-k, kh, kh, 0 
ies 10 12 13 21 Vy, 31 V, p 
V, ay, 
C. zelt = : 
where x = e , A= ky, fay 0 ; ,B=| 9 Jandc=[0 0 0 1] 
3 V, 0 
e hse 0 —k;, 0 ‘ 
3 
L ko 0 0 ~ Keg 





Ci, C2, C3 and C. denote the drug concentration in compartment 1, 2, 3 and effect compartment, 
respectively. V; is the volume for compartment i while 4; is the drug transfer rate from compartment i to j. u is 
the propofol infusion rate (ml/h), p is the drug concentration (for propofol, p = 10mg/ml) and a is a 
normalisation constant (a = 60 min/h). 

For Schnider PK model, the nominal value of compartment volumes and drug transfer rates can be 
found in [14-15]. However, due to inter-patient variability, it is expected that the true values will be different 
from the nominal values. Luckily, in most of the biological modellings, inter-variability remains bounded with 
a priori known bounds [16]. Consider a 50 year old male patient weighting 70kg and 170cm tall. A possible 
lower and upper bound of matrix A were given as the follow: 


—0.9689 0.3190 0.1958 0 —0.7979 0.3190 0.1958 0 
fz 0.0679 —0.0679 0 0 re 0.0679 —0.0679 0 0 o 
> 0.0035 0 — 0.0035 0 0.0035 0 — 0.0035 0 

0.2420 0 0 — 0.4820 0.2420 0 0 — 0.2420 


B and C were assumed certain and were equal to B = [0.039 0 0 OJ andC=[0 0 0 1|PD model can 
be expressed by the following static nonlinear Hill equation: 
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C’ 


BIS = BIS, +(BlS ya, — BIS, ) 
Clso + CZ 





(3) 


where BIS is measurable parameter of anaesthetic level, BISo is the base line effect in the absence of drug, 
BlISmax is the maximum effect of the drug, Ceso is the concentration that produces 50% of the maximum drug 
effect and y determine the steepness of the Hill equation. In this paper, it is assumed that BISo = 100, 
BISmax = 0, Ceso = 2.23 and y = 1.72. Due to the static nature of PD model, the nonlinear effect can be 
cancelled out by its inverse, PD"!. Hence, the plant can be controlled using linear control strategy. 


b. State feedback stabilisation 

First, it is assumed that the states are available for feedback. The objective is to determine the gain 
matrix K in the control law u = Kx that stabilizes the system = Ax+ Bu while guaranteeing the system 
positiveness. According to [8], the K can be determined such that A + BK is both a Metzler and a Hurwitz 
matrix. These can be achieved by the following LP inequalities: 


Ad+ BY. v, <0 
i=l 
Ad + BY'v, <0 (4) 
i=l 
-d<0 
And 
—v, <0 
- os vy, <0 
a (5) 
-a,d,—-byv, <0 fori j 
-a,d,-byv, <0 fori#j 
where K was computed as 
K =|d;"v, K d;'v, |. (6) 
c. Positive interval observer 
Consider the following observer system 
B= Ax 
y=Cx 7) 


where the trajectory of x € R” is assumed to be unknown but nonnegative. Based on the classical Luenberger 
observer [17], the following linear observer was used to approximate the states x 


&=(A-LC)x+Ly (8) 
where L is the observer gain. The aim here is to determine the observer gain Z such that the approximated 


states x are nonnegative and the error x —x converges asymptotically to zero. 
According to [10], a positive interval observer can be designed using the following inequalities: 


Indonesian J Elec Eng & Comp Sci, Vol. 10, No. 3, June 2018 : 1000 — 1006 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 O 1003 





ATd-C'Yz, <0 


i=l 


-d<0 
—c7z, <0 fori, 7=1, Kn (9) 
—A' diag(d)+C"Z-fI <0 
—z, <0 
with L calculated as 
L=|aj'z, K diz, |. (10) 


d. Integral action 
To include reference tracking and improve its disturbance rejection abilities, the control law was 
expanded to include integral action 


u=—Kx-k,z iy 
B=y-r (1) 


where r is the reference point and 4; as the integral term. The augmented system of the closed-loop controller 
was expressed as follows: 


A 0. OJ]x] [B 0 
=|LC A-LC 0/2|+/Blu+] 0 (12) 
C 0 Ofz}| jo} f-z 


Figure 1 shows the overall block diagram of the closed-loop system. The plant model has the 
propofol infusion rate, u as the input and Bispectral Index (BIS) as the output. BIS is the was assumed 
unknown. BIS is a widely used index that measures the depth of anaesthesia. The nonlinear effect introduced 
Since the patients were assumed unknown, the nominal matrix A was used in the observer block during 
simulation. 
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Figure 1. Block diagram of the observer-based state feedback control with integral action 


3. RESULTS AND ANALYSIS 
Based on the method described in Section 2, the controller gains were designed and computed to be 
K =[0.3593 0.3660 0.3759 0.5103], L=[0 0 0 0.0464], and k; = 2. The designed controller was 


tested on 9 simulated patients with varying ki9 and ko values, as tabulated in Table 1. These two values were 
chosen due to their higher impact on the variability of BIS effect. 

The objective of the controller is to bring the BIS value down to 50, which is the recommended 
index during surgical procedures. The controller was assessed by its induction phase duration (ID), 
percentage of overshoot (OS), and integrated absolute error (IAE) [18]. The ID was defined as the time 
elapsed from the beginning of propofol administration until the BIS falls to below 60 for 30s [19]. A short ID 
is preferable because it reduces the time spent in the operating room. OS measures the degree of overshoot 
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and should be as low as possible to prevent drug overdose. IAE is the integration of absolute difference 
between the BIS set-point and the real BIS. 


Table 1. Test Population 





Patient no. kyo (min!) keg (min'!) 
1 0.283 0.242 
2 0.283 0.362 
3 0.283 0.482 
4 0.369 0.242 
5 (Nominal) 0.369 0.362 
6 0.369 0.482 
7 0.454 0.242 
8 0.454 0.362 
9 0.454 0.482 





Table 2 and Figure 2 shows the simulation result of the controller for all patients. Result shows that 
all BIS converges to the targeted value. The mean value of OS is 13.01%, with no BIS lower than 40. 
However, a long ID, with the mean of 14.78 min, was observed. This duration was longer to 3 min, which is 
considered too long for clinical practice and other research [20]. 


(a) 


BIS 








Table 2. Performance of the controller without initial bolus 











Patient no. ID (min) OS (%) IAE 
1 14.50 19.42 13629 

2 13.50 14.70 11321 

3 13.08 12.34 10564 

4 15.50 16.33 12225 

5 14.58 12.08 10458 

6 14.25 10.12 9726 

7 16.50 13.77 11211 

8 15.67 9.98 9704 

9 15.42 8.37 9021 
Mean 14.78 13.01 10873 
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Figure 2. Performance of the controller without initial bolus. (a) BIS response (b) Propofol infusion rate 


The induction phase can be greatly reduced when an initial bolus of propofol was administered. Due 
to the presence of closed-loop control, a lower-than-clinical practice bolus injection was administered during 
the initial phase, i.e., 250 ml/h at the first minute. Then, the controller was initiated in the second minutes. 
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Table 3 and Figure 3 presents the simulation result of the controller with an initial bolus for all 
patients. The result shows that all BIS converges to the targeted value. The mean ID was reduced to 2.36 min. 


All BIS reading reached 60 within 3.03 min. However, the mean OS was shown to increase from 13.01% to 
16.60%. 


Table 3. Performance of the controller with initial bolus 





Patient no. ID (min) OS (%) IAE 
1 2.83 19.22 14047 
2 2.08 20.07 12855 

3 1.83 27.79 12291 

4 3.00 14.54 1843 
5 17 12.46 10520 

6 1.83 21.28 9944 

7 3.33 11.55 10477 

8 2.25 7.36 7975 

9 1.92 15.15 7871 
Mean 2.36 16.60 10869 
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Figure 3. Performance of the controller with initial bolus. (a) BIS response (b) Propofol infusion rate 


While the ID can also be shortened by increasing the 4; value, it will also cause a large overshoot 
which may put the patient safety at risk. Besides these two conflicting performance criteria (ID and OS), the 
performance of the controller has to be compromised for system robustness if the BIS response were to be 
kept within the range of 40-60 for all patients. This shows that the inter-individual variability is a great 
challenge in the design of a safe and efficient controller for anaesthesia system. 


4. CONCLUSION 

An observer-based state feedback control with integral action was designed to regulate the BIS 
signal using propofol infusion rate. The state feedback gain and observer gain were designed using the LP 
approach and took into consideration the positiveness of states as well as the uncertainty of the PK model. 
Simulation result shows that all BIS was able to converge to the set-point, but with a long induction phase. 
A bolus injection of propofol during the initial stage was able to reduce the induction phase, but will increase 
the overshoot percentage. This shows that the design of a fixed controller for the whole population is difficult 
due to the large inter-individual variability exists among patients. Individualised controller may be the key to 
tackle patient variability in order to realise the closed-loop control of anaesthesia. 
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